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Abstract 



^ . We compute tunneling in a quantum field theory in 1 + 1 dimensions for 

a field potential U(<&) of the asymmetric double well type. The system is 
localized initially in the "false vacuum" . We consider the case of a compact 
space (Si) and study global tunneling. The process is studied in real-time 
simulations. The computation is based on the time-dependent Hartree-Fock 
variational principle with a product ansatz for the wave functions of the 
various normal modes. While the wave functions of the nonzero momentum 
modes are treated within the Gaussian approximation, the wave function of 
the zero mode that tunnels between the two wells is not restricted to be 
Gaussian, but evolves according to a standard Schrodinger equation. We 
find that in general tunneling occurs in a resonant way, the resonances being 
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associated with degeneracies of the approximate levels in the two separated 
wells. If the nonzero momentum modes of the quantum field are excited only 
weakly, the phenomena resemble those of quantum mechanics with the wave 
function of the zero mode oscillating forth and back between the wells. If 
the nonzero momentum modes are excited efficiently, they react back onto 
the zero mode causing an effective dissipation. In some region of parameter 
space this back-reaction causes the potential barrier to disappear temporarily 
or definitely, the tunneling towards the "true vacuum" is then replaced by a 
sliding of the wave function. 
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1 Introduction 



Tunneling is one of the important elementary processes that may happen in a 
quantum mechanical system. There is a vast literature on the subject and the 
WKB approach is discussed in textbooks on quantum mechanics. In quantum 
field theory tunneling has been mostly discussed as a local process proceeding 
via bounce or bubble solutions P] 121 El H] for the simple reason that in an 
infinite space a global transition of a mean field through a barrier would 
have an infinite action. However, if space is finite such global transitions are 
possible. 

Tunneling with compact spaces can be of interest in cosmology. Such 
transitions in de Sitter space have been evoked in order to describe the quan- 
tum creation of a universe Ej. In a less specific way one may think of 
transitions within a finite volume of a chaotic initial state of the universe [7] . 
Tunneling in a finite volume may also have some relation to local tunneling 
in an infinite volume as an alternative to bounce transitions. This point has 
been elaborated in Ref. jH]. 

Tunneling with compact space also occurs in models for the quantum 
creation of a universe EI] , which in the minisuperspace space approxima- 
tion is a closed 3 sphere and whose radius a is a simple quantum mechanical 
degree of freedom. This effect has been widely discussed, and most authors 
use the WKB approximation ^T] Ell El CHI QUI- Our analysis is not 
directly related to quantum cosmology because the role of time is quite dif- 
ferent there and here. However, there are similarities concerning the role 
of particle production and their back-reaction which still are an open issue 
there|T71 ubi US HH] ■ 

Tunneling is of course a process that occurs in real time and the contin- 
uation to imaginary time is a technical tool which is widely used and whose 
application is based in a controlled way on the eikonal expansion. Its applica- 
tion to quantum field theory becomes cumbersome whenever one is interested 
in the physical evolution of the system during and after the tunneling pro- 
cess, the matching of the wave functions that is well understood in quantum 
mechanics now encompasses an infinite system of modes and becomes rather 
involved. It therefore seems to be of interest to look for a description of the 
system entirely in real time. 

Such an approach was taken for the case of quantum mechanics in Ref. 
[2Tij . These authors have studied in detail the case of the double- well poten- 
tial. They find that tunneling is characterized by the occurrence of resonances 
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between degenerate (approximate) levels in the spectrum of the separate 
wells. 

For quantum field theory the real time approach to quantum tunneling, 
in the same way as it is used here, has recently been discussed by Hirota [S] . 
His analysis, as ours, is based on the time-dependent Hartree-Fock approxi- 
mation. He uses analytic approximations for the zero mode wave functions. 
While quantum back reaction and renormalization are mentioned, they are 
not considered in detail and the author does not present any numerical re- 
sults. 

We do not solve the field theoretical problem exactly but use two ap- 
proximations: for the wave function of the entire system we make a product 
ansatz, i.e., a Hartree approximation. The wave functions of the fluctuation 
modes are taken to be Gaussian; for the zero mode wave function, however, 
we numerically solve the Schrodinger equation. In this respect we differ from 
the widely used out-of-equilibrium simulations [2U 122 I2H1 EH EH] , where a 
Gaussian wave function, shifted by a classical field </>(£), is used for the zero 
mode as well. By using the Hartree approximation we go beyond the one- 
loop back-reaction by taking into account the back-reaction of the nonzero 
momentum modes onto the zero mode and onto themselves. The latter back- 
reaction avoids a possible instability of the system whenever the effective 
mass of the nonzero momentum modes gets imaginary. 

The plan of the paper is as follows: in Sec. [21 we present the model, the 
decomposition into a discrete set of degrees of freedom, and the Hamiltonian 
in the Schrodinger representation; in Sec. |3]we formulate the time-dependent 
Hartree approximation [2H1 123 I2H] using our ansatz for the product wave 
function; renormalization is discussed in Sec. E]and the initial conditions in 
Sec. EJ numerical results for some parameter sets are presented and discussed 
in Sec. El we end with conclusions and an outlook in Sec. [7| 

2 The model 

We consider a scalar quantum field theory in 1 + 1 dimensions with a com- 
pact space, i.e., on a manifold R x Si. As ultimately we want to use our 
approach on R x S% we only introduce those self-couplings that would lead 
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to a renormalizable theory in 3 + 1 dimensions. The action is given by 



2na 



dx / dt 



1 /<9$y 

2 Vary 



with 



[/($) = i m 2 $ 2 - ??$ 3 + -$ 4 
2 8 



(2.1) 



(2.2) 



and periodic boundary conditions in x. With positive values of the couplings 
the potential has the asymmetric double well form with a minimum at $ = 
$_ = and a second one at 



1 a 
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(2.3) 



with a = \m 2 /4r) 2 . The parameter a is restricted to the range < a < 1. 
The field $ is dimensionless, and rj and A have the dimension mass 2 . We 
introduce the expansion into normal modes 



$(x,t) = (p (t) + V"(*) e 



(2.4) 



with 



Then the action takes the form 



n 



S = 2na I dt -(f 2 Q + ^ ^ ¥>n¥>-n ~ ^ J^^-n 



with 



A x - 

g 2^1 ( Pri ( Pn> ( Pn" l P-n-n'-n" 
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or 



(2.5) 



(2.6) 



(2.7) 



As the field $ is real we have the condition (p_ n = <£>t . We therefore introduce, 
for n 7^ real fields via 



f±n = ~^P'^ ± i( Pn2), n > . 



(2i 



4 



In the Schrodinger representation the Hamiltonian is given by 

1 



H 



2na 
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+ 2 
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of V — 
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m (f 



5>v 



nrnj 



n>0,j 



-V ^ ^n^n'^-n-n' + ^ ^ y?nWW'y? -n-n> -n" 



(2.9) 



In order to get rid of factors 27ra we introduce the rescaling 



(f k = Xfc/v27ra , 
so that the Hamiltonian takes the form 



H 



1 

+ 2 



d 2 







2 2 X ^ 2 2 



n>0J 



A' 



^7 ^ ^ XnXn'X—n—n' ~l~ q ^ ^ XnXn'Xn"X— n— n'— n" 



(2.10) 



(2.11) 



with Vj — r\l and A' = \j2ixa. Note that in the cubic and quartic parts 
we still have retained the complex fields 



X±n = ■J=(Xnl±iXra) ■ 



(2.12) 



3 The time dependent Hartree-Fock approx- 
imation 

We now make a variational ansatz for the wave function and impose a varia- 
tional principle, which is known under the name of "time-dependent Hartree- 
Fock approach". The ansatz for the wave function is 



*(X0, Xn, t) = 1p (X0, t) Y[ ^n(Xnj, t) 



(3.1) 



n>0 
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Furthermore we will restrict the ansatz for the modes with n to a Gaus- 
sian wave function 



-ia n (t) 



[2*o*(t)]V* 



exp 
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— i- 



■ S n (t) 



2 \2al{t) a n (t))^ 



A 7. 



(3.2) 



while we do not further specify ipQ. The time dependent variational principle 
[231211], now imposes the condition 



5 MXo [J dxnj^(Xo,Xnj,t) (id t - H) %(xo,Xnj,t) = . 

n>0,j 



(3.3) 



We will find later that the dynamics is independent of j, which therefore has 
already been suppressed in the index for the wave functions. Before we write 
down the resulting Schrodinger equations for the various degrees of freedom 
we rewrite the Hamiltonian in an appropriate way. 

The part of H which exclusively contains the zero mode is given by 

1 d 2 

Hm = -^— + U{xv) , (3.4) 



with 



U( X o) = \m 2 xl - ixl + jxt • (3-5) 
The parts bilinear in the quantum modes lead to the Hamiltonian 

1 d 2 1 2 2 " 

2 8x1, ^ nXn \ 



On 



E 

n>Q,j 



(3.6) 



Due to the Gaussian ansatz for the n ^ modes the expectation values of 
any odd powers of n ^ fluctuations vanish. So if we consider the interaction 
between zero and nonzero modes we need to retain only terms bilinear in the 
n ^ modes. We then have for the interaction between zero and nonzero 
modes 



/ , XnXn'X~n—n' 



3X0 XnX-n = 3xo 



X 



"j 



(3.7) 



70n 



n^O 



n>0,j 



and 



E 



XnXn'Xn"X—n—n'—n" 



6Xo XnX ~ n = 6 *o Yl X*i ■ ( 3 ' 8 ) 



/On 



n^O 



n>0,j 
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So with the scaled couplings rf and A' the interaction Hamiltonian cou- 
pling zero and nonzero modes is given, in the approximation considered here, 
by 



3A' 

H IOn = (-3?/xo + ~^-xl) Yl 



X 



nj 



(3.9) 
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The part of the Hamiltonian cubic in the nonzero modes has a vanishing 
expectation value. The quartic term has the expectation value 



< Yl XnXn'Xn"X-n-n'-n" >= 3 ^ < xlj >< xl'j' > , 



(3.10) 



where we have used the fact that the wave functions for n^O are Gaussian. 
This part yields the Hamiltonian for the mutual and self interaction of the 
nonzero modes. It is given by 



< H Inn >— -A' 



nj 



(3.11) 



Applying now the variation principle we get for ip the Schrodinger equation 



id t ipo(xo,t) 



Id 2 ~ 3A' 
2 + U(xo) + (-6t/xo + — Xo)-^ 



Here we have introduced the fluctuation integral 



^ = 2 < Y Xlj > , 

n>0,j 



Mxo,t). (3.12) 



(3.13) 



which will be specified more explicitly later. For the modes with n^Owe 
find 



idti> n (Xn,t) 



1 d 



2 +-(u; 2 n + W(t))x 2 nj 



Here 
with 



W(t) =< U"( X o) -m 2 > +3A'J^ , 

3A' 



ll>n(Xn,t) 



(3.14) 

(3.15) 



< U"( X o) -m 2 >= / d X or o (Xo,t)[-6v'Xo + -^xl}Mxo,t) . (3.16) 



We also introduce fi^(t) = u> 2 +W(t). The Schrodinger equation for the quan- 
tum mode wave function is seen to be independent of the index j. Therefore 
the previous summations over j just lead to a degeneracy factor 2, so that 
for the mode functions the index j can be suppressed. 

With the Gaussian ansatz ()3.2|) the Schrodinger equation for the n ^ 
modes implies 



Sn(t) 

a n (t) 



Sn(t) , 

1 



(3.17) 
(3.18) 

(3.19) 



The first two of these equations can be related to mode functions f n (t) sat- 
isfying 

f n (t) + n 2 n (t)f n (t) = 0, (3.20) 

as they arise from a Klein-Gordon equation for a field with the effective mass 
m 2 + W(t). We have, with u n0 = fi n (0), 

!/»(*)! 



<?n{t) 
S n (t) 

while the wave function is given by 



1 d 



\J1oj n Q dt 



\fn(t)\ 



(3.21) 
(3.22) 



-ia n (t) 



2u. 



1/4 



exp 



(3.23) 



In deriving these relation one uses repeatedly the Wronskian relation 

£(*)/«(*) - £(*)/»(*) = 2iu n0 , (3.24) 
which corresponds to an initial condition 

/ n (0) = 1 , / n (0) = -iu n0 . (3.25) 
For this wave function the expectation value of x\ * s given by 

\fn(t)\ 2 



2 ,,2 
n 



dXn\*Pn(Xmt)\ X 
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2u. 



nO 



(3.26) 



so that the fluctuation integral becomes 



WW 



n>0 



2u. 



nO 



(3.27) 



The system of equations we have presented here is consistent with an 
conserved energy 



E =< H o0 + H 0n + H IOn + H Inn > 



with 



< #00 > 

< Hon > 

< Hio n > 

< Hi nn > 



1 d 2 



+ U{Xo) 



6-*n 2a; n0 



n>0 

< U"( X o) -m 2 >T 
3 



fn(t)f + lu 2 n \f n (t)t 



4>o(Xo,t) 



/-r-2 



-XT 



(3.28) 

(3.29) 
(3.30) 

(3.31) 
(3.32) 



It is convenient to replace u 2 by u 2 in < H 0n >, so that < H 0n > becomes 
the free Hamiltonian for the initial fluctuation wave functions. Then the 
term < H IOn > receives an additional contribution — iy(0)jF. 

A particle number for the n ^ modes may be defined in various ways. 
A Fock space is defined by the mode decomposition ()2.4|) . To the operators 
Xnj and their conjugate momenta n n j = d/dxnj we can associate creation 
and annihilation operators referring to an oscillator of frequency u> n o via 



- "J 



UnO ( 1 9 

2 \ Xn:> U n0 dXnj 

1 d 



^>n0 



Xnj 



^nO 9Xnj 



and a particle number 



"J 



^ C nj C nj > — < 
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1 d 2 1 2 2 

,2 i ry^noXnj 



2d X 



> 



(3.33) 
(3.34) 

(3.35) 



computing the expectation value with the wave function ip n (Xnj), Eq- ()3.23|) . 
one finds 



n.) 



f n (t) +U 2 n0 \f n (t)t 
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1 
2 



(3.36) 



In the present context the quantum excitations do not describe free particles; 
we merely use the total "particle number" J\f = J2 n j N n j as an indicator for 
the excitations of the n > oscillators. 

4 Renormalization 

The fluctuation integral and the energy density, as introduced in the previous 
section, are divergent quantities. In 1 + 1 dimensions the renormalization is 
in principle rather straightforward. However, here it has to be done for a 
nonequilibrium system with discrete momenta, and for the Hartree approxi- 
mation, i.e., for a resummed perturbation series. We here limit ourselves to 
present the main steps needed to derive the formulas that enter the numerical 
codes. 

Though we treat the system in a nonperturbative approach, the diver- 
gences are related exactly to those of standard perturbation theory. The 
mode functions can be expanded perturbatively with respect to the poten- 
tial V(t) which contains the couplings 77 and A; so such an expansion is at 
the same time an expansion with respect to these couplings. 

As discussed in previous publications on nonequilibrium dynamics ;24 
we write the mode functions as 



and convert the differential equation for the f n (t) into an integral equation 
for the functions h n {t) as 



Ut) = e 



[l + h n {t)] 





) - l) V(0 [1 + K(t')) . 



(4.2) 



Here V(t) = W(t) - W(0), see Eq. (j3~T51t . 



We further have 



\f n (t)\ 2 = l + 2Reh n (t) + \h n (t)\ 2 . 



(4.3) 



One easily finds that Keh n (t) behaves as 




(4.4) 
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for large n. For the fluctuation integral one finds 



^ 2u; n0 



2^o 



10, 



iM 



(4.5) 



so that we can separate it into a divergent and a subtracted finite part as 

n=0 
oo j 

-^sub = J2o [2Re/i n (t) + |/i n (t)| 2 ] . (4.7) 



n=0 



nO 



The divergent part has to be regularized and to be separated into the stan- 
dard renormalization part and a finite contribution. The sum over discrete 
momenta here and below can be done in the same way as the Matsubara 
sums in finite temperature quantum field theory, or using the Plana formula 
[HI. One finds 



oo 1 

V — 



na 



dk 



_ 00 (2vr)2u;o 



1 + 



exp(27rac<Jo) — 1 iracvo 



(4.8) 



with oo = \Jk 2 + m\. The first term in the bracket is obviously the divergent 
part we would obtain for infinite space and goes into the renormalization of 
the various couplings. The second term in the bracket arises from the periodic 
boundary conditions. The third part arises from the fact that T contains the 
nonzero modes only and no subtraction has been applied to the zero mode. 
The first term can be regularized in a standard way [21]; we have 



dk 



(2tt)2o;c 



rex 



d 2 ~ f k 



(2tt) 2 ^ kl - k 2 - m 2 



+ io 



1 r 

4^°' 



(4.9) 



with L t Q = 2/e+ln(47r/mQ)— 7 using dimensional regularization. m\ depends 
on the initial conditions, and the renormalization condition should not; we 
therefore write L e q = L e + \wm 2 jm\ with L e = 2/e + ln(47r/m 2 ) — 7 and 
include the second term into the finite part. So we define the remaining 
finite part as 



(o) 



fin 



2ixa 



dk 



1 



(2tt)2oj exp(27raa;o) 



a m 
4 mi 



(4.10) 
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and the finite part of the fluctuation integral as 



(0) 



sub 



We similarly decompose the fluctuation energy as 

'1 A , 2 1 



< H(\ n > 



2 > — 



n>0 



1 



.un + ~u 2 n \fn(t)\' 



(4.11) 



(4.12) 



-Efi.sub - -m + a I dkuo - 2a 



k 2 dk 



uiq (exp(27raco> ) — 1' 



with 



fl.sub 



2E — 



n>0 



/n(*) 



| w nl/n(*)| 2 - 



(4.13) 



Integrating by part the last ( "thermal" ) integral in Eq. (|4.12jl can be recast 
into the form of a free energy 



-2a 



k 2 dk 



dk 
27 



ln(l - exp(-27rac«Jo)) . (4.14) 



uiq (exp(27rau; ) — V 
The integral over uoq is given, in dimensional regularization, as 

,2 



dkuo 



o 



r<'S 



m 



-m n [ L e + In — o + 1 



(4.15) 



and we define 

Efi M = £ fl , sub + -m 2 , In — + 1 + / -In 1- exp -2vraw . 4.16 
4 V. mo J Jo 2k 

The renormalization is done in analogy to the case of nonequilibrium 
dynamics in , following the scheme of Ref . j22] , by adding a counterterm 

SCM 2 = 6C(m 2 + < U"( X o) ~m 2 > +3\'^ n ) . (4.17) 
Here one has used already the finite gap equation 

M 2 = m 2 + < U"(xo) ~m 2 > +3\'F Rn , (4.18) 

which determines the effective mass of the fluctuations. Choosing 

5C = ~L e , (4.19) 
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and denning the initial mass m by the initial gap equation 

ml = m 2 + < U"( Xo ) - m 2 > +3X'J^ , (4-20) 

the gap equation (|4. 18|) is satisfied for all t. It can be written as 

M\t) =m 2 + V(t) , (4.21) 

with the previously introduced potential 

V(t) =< U"(xo) > -m 2 + 3A% n (t) , (4.22) 

which now is well-defined. 
For the energy we obtain 

1 3 

E = < Foo(xo) > +E RM - -m + < U" - mg > + -A'J^ s 2 ub 

a m 2 f 00 dk 

+-m 2 (ln^ + l)+ / — ln(l -exp(-2vra^o)) . (4.23) 
4 ml J 2tt 

The last two terms are, for a fixed radius a, independent of time; we leave 
them out when presenting the energy conservation. The term —mo/2 marks 
the absence of the zero mode in the sum over fluctuations. If mo = m this 
term exactly cancels, at t — 0, the expectation value < if o(Xo) >• Likewise, 
the term — l/4mo i n the initial gap equation ()4.20j) is cancelled, for mo = m, 
by the term < U"(xo) ~ m ^ >• Indeed if a is sufficiently large, am > 1, 
corresponding to a "temperature" T e g/m < 1/27T, the "thermal integrals" 
become very small; then the solution m of the gap equation indeed is close 
to m and E w 0. 

5 Initial conditions 

As already discussed in Sec. OJthe system is started for the nonzero modes 
with 

f n (t) ~ e—* , (5.1) 
which is equivalent to initial wave functions 

e -^noxl/2 . (5.2) 



7T 
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Here u n0 = yrr^ + n 2 / a? and itlq is determined by the gap equation (j4.2()j) . 
For the zero mode we likewise start with a Gaussian wave function 



V>(xo,o) 



mi V4 



e -m X l/2 _ 



(5.3) 



This would be the ground state wave function of an oscillator with frequency 
uj = m, i.e., the ground state wave function for a potential with 77 = A = 0. 

The choice of these initial conditions is of course quite arbitrary. Here 
we try to start in an approximate ground state for the system "in the left 
well" of the double well potential. We find the total energy of the system to 
be close to zero, as one would expect for such an approximate ground state 
(it is understood that the zero point energies of the nonzero modes are not 
included). A different choice would in general lead to a higher total energy 
and make tunneling easier, or make the transition to a large part an "over 
the barrier" transition. 

6 Numerical results 
6.1 Numerical details 

We have implemented the formulas of the previous sections numerically. This 
is essentially straightforward. The Schrodinger equation for the zero mode 
becomes a system of first order differential equations for the values "0o(Xo,fc) 
where Xo,k are the equidistant discrete values of the variable \o, an d the 
coupling arises from the discretized second derivative 



where Axo is the step width for the xo,k- The Schrodinger equation for the 
n 7^ modes is converted into the second order mode equations, which are 
coupled to the zero mode and the other nonzero modes by V(t). This dis- 
cretization leads to instabilities in the time evolution unless the time intervals 
are chosen of the order A%q. We used 4000 grid points for the zero mode 
Xo, which typically extends over a region of —10 < xo < 50. So Axo is of 
the order 10 -2 . The Runge-Kutta time step was chosen At = .00002. This 
choice makes the codes very slow, much slower than those of nonequilibrium 
dynamics in the Gaussian approximation. As we do not use the Gaussian 



^o(Xo,fc) 



ip(Xo,k+i) + ip(xo,k-i) ~ 2^(xo,fc) 
(Axo) 2 



(6.1) 
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Figure 1: Energy conservation; parameters m — rj — X — 1, a — 6; 
dotted line: classical energy, Eq. (j3.29j) : dash-dotted line: fluctuation en- 
ergy, Eq. (j3.30ft : dashed line: interaction energy, Eq. (j3.3H) : solid line: self- 
interaction of the fluctuations, Eq. ()3.32|) . The total energy (not displayed) 
is E = (2.39 ±0.001) x 10" 5 . 

approximation for the zero mode wave function, we have to compute, at each 
Runge-Kutta step, not only the sums over the quantum modes \m Du t also 
the averages of various observables in the ground state wave function. The 
excited quantum modes were taken into account up to n = 200. A simulation 
of the time evolution until t ~ 100 takes a few hours on a standard PC. 

The accuracy of the computations was checked by computing Wronskians 
and energy conservation. The relative accuracy obtained was better than six 
significant digits. The energy conservation also checks the correct implemen- 
tation of the basic equations. We show a typical example, for A = r\ — 1 
(set IV below) and a = 1.2. The total energy is E = (2.39 ± 0.001) x 10~ 5 
throughout the total time interval, the single components take values up to 
35. 
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6.2 Parametrization and parameter sets 

The parameter space, encompassing 77, A and a is quite large, so before pre- 
senting numerical results we try to get some qualitative insight into the 
physics to be expected in certain ranges of their numerical values. The pa- 
rameter m sets the overall scales and is put equal to m — 1. 

While in Sec. |2l we introduced a scaling of the fields, (|2.1(J|) . which was 
suitable for the canonical formalism on the basis of the Hamiltonian ()2.11|) . 
the discussion of the results is more transparent if we introduce the a — (3— 
parametrization widely used in bounce computations |U |33J |3U One 
introduces the rescaling X = mx and $ = m 2 $/2?7. Then for infinite space 
the classical action takes the form 



s cl = ^ 2 Jd 2 x [I (v,$) 2 + f >($) 



(3S cl (a) , (6.2) 



with (3 = m^/Ar] 2 , a = X/3/m 2 and 



U(§) = i$ 2 - H 3 + -<f> 4 . (6.3) 
2 2 8 

While f3 multiplies the classical action the one-loop effective action is a func- 
tion of a only. So for large (3 the system essentially becomes classical and the 
quantum effects only lead to small corrections while for small (3 the quantum 
effects become large. The parameter a determines the shape of the potential: 
for a —>■ 1 we get a symmetric double-well potential, for small a the right 
hand well becomes much deeper than the left hand one, the barrier between 
the wells becomes shallow. For bubble nucleation a — * 1 is the thin-wall 
limit, where the bubble size becomes very large. The case a > 1 is of no 
interest here, as the minimum at $ = then becomes the global minimum. 

We will in the following consider parameter sets with fixed a and /3, i.e., 
fixed A and rj, and study the dependence on a. 

For finite space extension the semiclassical tunneling action of a spatially 
homogeneous bounce is given by 

bounce = 2rraP [ ° d$J&(&) , (6.4) 

where 



$0 = - (1 - VT^) (6.5) 

a 
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is the zero at the right hand side of the potential barrier. For fixed a, (3 the 
tunneling via a homogeneous bounce should shut off with increasing a. The 
tunneling will then occur via local bounces, as considered recently in the 
Hartree approximation in |37)1 I3~T] . 

The transition rate obtained from the homogeneous bounce, Eq. ([6.4)1 . 
will only have a qualitative relation to the observed tunneling phenomena 
which will be characterized by the occurrence of resonances. For quantum 
mechanical tunneling this was already observed in Ref.j2HI- I n order to esti- 
mate the separation of the resonances we consider the approximate spectra of 
the left and right wells, treating them as separate oscillators. The resonances 
can then be thought of as arising from the degeneracy of levels in the left 
and in the right wells. From the Hamiltonian and the potential written in 
canonical variables, Eqs. (|3.4|) and (|3.5jl . we see that the energies of the left 
hand oscillator are E l n = (n+ l/2)m. The right hand oscillator potential has 
its minimum at 



with 



Xo 



A 




(6.6) 



(6.7) 



and the energy levels are given by 



El 



(6.8) 



Here 



and 



<4 = U"(xt) = 2m 2 




(6.9) 



U{xt) = 2na(3U(^ + ) < . (6.10) 

So the condition for a degeneracy of a level in the spectrum of the right hand 
well with the ground state of the left hand well (our initial state) is 



2vra[/($ + ) + ( n + - ) oo A 



1 

—m . 

2 



(6.11) 



We present here the dependence of the tunneling phenomena on the spatial 
length scale a. If the degeneracy holds for some value a and an integer n, it 
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will hold again for a' = a + Aa, n' = n — An with 



Aa = 



An . 



(6.12) 



2nU(<f>+) 



The constant multiplying An on the right hand side determines the spacing 
of resonance levels as a function of a at fixed a and (3. Of course we may also 
have resonances between excited states of the left well and those of the right 
hand well, but, as our initial state will roughly correspond to the ground state 
of the left hand well, these will be less important. Another, more essential 
feature is the excitation of field quanta of the nonzero modes. These will 
have a dissipative effect on the dynamics of the zero mode and broaden the 
resonances. The interaction consists in a deformation of the potential in 
which the zero mode is moving, which takes the form 



If J- it) is negligible the evolution of the system proceeds like in the quantum 
mechanics of the zero mode. If T(t) remains small this time-dependent 
modulation of the potential will allow a few other approximate eigenstates of 
the zero mode to mix in, the resonant oscillations develop "higher harmonics" 
and become irregular. If jF(t) gets large then the potential can be deformed 
in such a way that the potential barrier disappears entirely, in such cases 
the zero mode may "slide" into the new minimum. This happens if J-{t) is 
positive; if F{t) is negative the potential is tilted counterclockwise and the 
barrier is enhanced. 

6.3 Results of the numerical simulations 

We have performed a study of tunneling as a function of the radius of the 
space manifold Si, for fixed values of the parameters a and (3, or A and n; 
the mass is chosen to be unity, which determines length and time scales. 
We have considered four parameter sets: set I: a = 0.8,(3 = 0.5, i.e., A = 
1.6, n = l/y/2; set II: a = 0.6,(3 = 2, i.e., A = 0.3, 77 = 1/V8; set III: 
a = 0.4, (3 = 1, i.e., A = 0.4, 77 = 1/2; set IV: a = 0.25, (3 = 0.25, i.e., 
A = 77 = 1. For set II we expect the quantum corrections to be small, for 
sets / and IV we expect large quantum corrections, and moderate ones for 



We have studied in general the behavior of the average of the zero mode 
< Xo(t) >; °f the fluctuations (fluctuation integrals, particle number, and 



U (xo,t) = U( X o) + (U"(xo) ~ m 2 )F{t) ■ 



(6.13) 



set m. 
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the various energies), and of the wave function ipo(xo) °f the zero mode, as 
functions of time. As the main indicator of the general behavior we use the 
maximal value attained by the expectation value of the zero mode during the 
time evolution. We denote this value as (po- In those cases where we have 
effective tunneling, this zero mode average settles at a value of ipo beyond 
the potential barrier, the late time average being typically 20% smaller than 
the maximal value. 

We display in Figs. 1 to 4 the dependence of the maximal values y?o of 
< tpo(t) >, as functions of a for the four parameter sets. The horizontal lines 
labelled by $ m and $o indicate the position of the maximum of U(<&) and the 
zero of the position between the two minima (the "end of the tunnel"). We 
see in all cases that an effective tunneling occurs as a resonance phenomenon. 
According to our estimate in the previous subsection, Eq. (j6.12|) . the spacing 
of different resonances as a function of a would be given by Aa = 0.529 for 
set I, Aa = 0.0338 for set II, Aa = 0.0167 for set III and Aa = 0.0166 
for set IV. The observed spacings roughly correspond to these estimates. 
We note that also the excited states of the left hand oscillator or of the 
nonzero modes can come into play, so if we observe some irregular spacings 
this may be due to such effects. When compared to similar figures obtained in 
quantum mechanical simulations [20] the resonances seen in our simulations 
are broader. We observe that even on resonance the maximal value of < (po > 
never reaches the second minimum (the "true vacuum"). In part this is due 
to the fact that the effective potential seen by the zero mode is modified by 
the quantum fluctuations, moreover, even at late times the wave function 
generally retains a finite probability density near the false vacuum tpo = 0. 
This will be discussed in more detail below. 

For the parameter sets / — III the tunneling shuts off entirely at higher 
values of a as expected from the bounce action (J6.4j) being proportional to a. 
As large (3 implies a large classical action and, therefore, a small semiclassical 
tunneling rate, this shutting-off already happens at small values, a ~ 0.5, for 
set II. For set / on still finds resonant tunneling for a > 3. However, 
the transition time at the resonances increases substantially (t ~ 1000 at 
the last peak in Fig. 1), between the resonances the amplitude of the zero 
mode oscillations decays exponentially with a, as for sets // and and 
the zero mode wave function displays no tunneling at all. The case of set 
IV is quite different. Here at large a the system always tunnels. This effect 
is due to the fact that the quantum fluctuations deform the potential seen 
by the zero mode, Eq. (J6.13|) into a potential without barrier, the late time 
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Figure 2: The maximum <p of the expectation value < ifo(t) > for set I, 
a = 0.8, = 0.5, function of a. 
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Figure 3: The maximum <p of the expectation value < <po(t) > for set //, 
a = 0.6, (3 = 2, cLS cL function of a. 
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wave function extending over a range from —5 to +20, we call this kind of 
transition a sliding transition. We will discuss this below. 

The results for the maximal value <po of the zero mode average display 
a global picture of the tunneling in the asymmetric double well potential 
in quantum field theory. They look quite different from what one usually 
assumes to happen in "false vacuum decay", when one just considers the 
homogeneous bounce action (J6.4)) without or with quantum corrections. We 
will now discuss in some more detail the evolution of the zero mode wave 
function. 

For all parameter sets and all values of a the wave function of the zero 
mode initially evolves slowly; it becomes slightly asymmetric and develops a 
tail into the region of the potential barrier. Here already the fact that we do 
not restrict it to be Gaussian is of prime importance. The further behavior 
then depends on the parameter sets. 

On the resonances the wave function, which initially is essentially the 
ground state wave function of the left well first connects to a particular 
wave function in the right well, obviously the one of the degenerate level. 
This is seen in Fig. El , the well-defined number of peaks within the right 
hand well indicate a specific approximate level within the right hand well. If 
the fluctuations of the quantum field, i.e., the amplitudes of the <p n , remain 
small, then one observes an oscillation forth and back between these two 
wave functions, as expected in a purely quantum mechanical system. If the 
nonzero field modes get excited, then these oscillations of the zero mode are 
disturbed, other wave functions mix in, as seen by the increasing number of 
peaks, in some cases the wave function looks quite chaotic. 

Off resonance the expectation value of <po oscillates regularly with a small 
amplitude. It still develops a tail into the right hand well, but this part of 
the wave functions remains small. 

Beyond the resonance region there is efficient tunneling if the fluctuations 
of the nonzero momentum modes are large. In this case, while the wave 
function penetrates into the barrier region, the fluctuation integral becomes 
positive, the barrier disappears and the wave function slides down the hill 
retaining almost its initial Gaussian form. 

The simulations of set II, a = 0.6,(3 = 2 are expected to resemble most 
closely the case of quantum mechanics. The nonzero momentum quantum 
fluctations indeed remain small and we have the resonances expected from 
the qualitative picture of degenerate levels of the individual wells. The re- 
sults for the expectation value of the zero mode and for the evolution of the 
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Figure 6: Probability density IV'oCxo) | 2 at intermediate times, here t = 10, 
for a = 0.8, (3 = 0.5, and a = 1.35 

wave functions correspond to these expectations. We present the regular os- 
cillations of the zero mode on the resonance at a = 0.5 in Fig. [7| We plot, in 
the same figure, the particle number Af defined below Eq. ()3.36|) . multiplied 
by 1000. The wave function, shown at four different times in Fig. |H1 is seen 
to return to itself almost exactly after half a period of t ~ 55. 

Set III is intermediate in the sense that with (3 = 1 we expect the fluctu- 
ations of the nonzero modes to be sizeable but not really large. We find that 
tunneling is again characterized by the occurrence of resonances, and again 
shuts off at larger values of a, here around a = 2. Around this value the 
resonant transitions occur on longer and longer time scales. So for a > 2 we 
cannot exclude further resonances if we run the simulations for times longer 
than t ~ 300. 

For set IV the zero mode shows resonant behaviour at small a, but for 
a > 0.5 the system always ends up in the right well. The actual evolution 
is quite involve here. Once the fluctuations set in, they tilt the potential 
in such a way that the barrier disappears and that the wave function can 
start to move right. Then the fluctuation integral T gets small again or 
negative, so the potential barrier appears again. This process repeats itself 
in an oscillatory way, and the wave function gradually shifts towards the 
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Figure 7: Evolution of the expectation value of the zero mode < (p > an d 
of the particle number Af; parameters a = 0.6, j3 — 2, a — 0.5. 
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Figure 8: The zero mode wave function for a = 0.6, (5 — 2, a — 0.5; dashed 
lines: the wave functions at t — and 109; solid lines: the wave functions at 
t = 50 and t = 168. 
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Figure 9: Evolution of the expectation value < (po(t) > an d of the fluctuation 
integral F{t) for a = (3 = 0.25 (set IV) and a = 1.2. 

deeper well. This is displayed in Fig. At late times the quantum modes 
still oscillate with a sizeable amplitude, so that the potential seen by the zero 
mode keeps on changing and can neither be considered to be double or single 
well. A typical wave function at late times is shown in Fig. El for a = 1.2 
and t = 30. One sees no trace of a double-well structure. In spite of the 
complexity of the wave function the energy is still conserved better than one 
part in 10 8 (see Fig. |TJ). 

The situation becomes more transparent for larger values of a 3 . We show 
results for a = 6. Here again the potential is tilted by the fluctuations, but 
the fluctuation integral remains positive on the average at late times. This 
implies that the barrier has disappeared definitively. The transition can be 
described as a sliding of the wave function towards the new minimum. This 
is displayed in Fig. [TT] 

For the case under consideration, set IV with a = 0.25 the potential is 
very asymmetric. So one may infer that the tilting of the potential is most 
effective here, and that this effect may in fact be limited to small values of 
a. We therefore have performed simulations for a = 0.8 and j3 = 0.2, i.e., for 
a more symmetric potential, but again with large quantum corrections. We 

3 These values do not appear in Fig. 0as otherwise the structure in the resonance region 
cannot not be resolved. 
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Figure 10: The wave function after tunneling, a = f3 = 0.25 (set IV) and 
a = 1.2 r~ A * - Qn 
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Figure 12: Evolution of the expectation value < <po(t) > and of the fluctua- 
tion integral for a = 0.8, (3 = 0.2, and a = 5. 

find that for values of a > 1 the transition again occurs by sliding, but that 
at late times the fluctuation integral becomes negative and the potential is 
tilted in such a way that the barrier reappears. This is presented in Figs. 
IT2*1 and The wave function at late times is then found to have sizeable 
amplitudes in both wells. 

So it seems that sliding occurs universally for small (3 and large a. 

7 Conclusions and outlook 

We have presented here an analysis of global tunneling in quantum field 
theory on a compact space. This analysis was based on a real-time formula- 
tion, using the time- dependent Hartree approximation. We have performed 
numerical simulations of a system where the wave function of the zero mo- 
mentum mode evolves as a solution of the Schrodinger equation, while the 
nonzero momentum modes are treated in the Gaussian approximation. The 
approximation includes the back-reaction of the nonzero momentum modes 
onto the zero momentum mode and onto themselves. The renormalized equa- 
tions were obtained using the 2PPI formalism, adopted here to a system with 
finite space extension. 
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Figure 13: Evolution of the effective potential U (xo,t)] solid line: t = 0; 
dashed line: t = 11; dash-doubledotted line: t = 13. 

We have found that tunneling for such systems occurs in a variety of 
different ways. For large (3, implying small field fluctations, the system be- 
haves like a quantum mechanical system with a single degree of freedom. 
Tunneling is effective whenever it connects degenerate modes, so varying the 
parameters, here the radius a, one finds resonance enhancements. For large 
spatial extension a this resonant tunneling shuts off as expected from the 
fact that the WKB integral over the barrier is proportional to a. For smaller 
values of (3 the nonzero momentum modes, the modes which make up the 
quantum field, are enhanced and modify the quantum mechanical behavior. 
If this enhancement is weak it can be considered as a kind of dissipation; the 
mean value of the zero mode shows regular periodic oscillations, on and off 
resonance. On resonance this again involves two or a few more modes in the 
two separate wells, off resonance the wave function essentially remains in the 
left well, with some tail in the right hand well. For parameter sets where the 
excitation of the nonzero momentum modes becomes sizeable, these regular 
oscillations are disturbed, the oscillations of the mean value of the zero mode 
become irregular and so do the wave functions of the zero mode. Off reso- 
nance at late times the system again remains concentrated in the left hand 



28 



well. On resonance at late times the wave function extends over the entire 
region allowed at the given energy. While the system may then be considered 
as having tunneled, it certainly ends up in a rather complex excited state. 

If the quantum corrections are large the system exhibits another phe- 
nomenon: the quantum fluctuations tilt the effective potential for the zero 
mode in such a way that the barrier disappears entirely. In such cases the 
system wave function can be observed to slide into the new minimum. The 
disappearance ot the barrier may be a temporary effect. It is due to the 
substantial quantum excitations once the wave function enters the barrier, 
leading to negative squared masses for the ffuctations. At late times the 
potential may become be a double-well potential again, or the barrier may 
disappear defintively. 

This variety of phenomena observed here in a real-time analysis cannot 
be expected to be described by a simple transition rate formula. A better 
approximation than the classical bounce action, Eq. ([6.4)1 . could be obtained 
by fitting together WKB patches within the allowed and forbidden regions 
or by applying a time-dependent WKB approximation to the k = mode 
itself, as done for quantum mechanics in Ref. |2Zj. However, the advantage 
of dealing with analytic or semi-analytic formulas is lost at the latest if one 
includes the nonzero momentum fluctuations. In most cases they cannot be 
evaluated analytically and, as implied by the present analysis, they cannot be 
expected to give reliable estimates. For homogeneous tunneling they exhibit 
the phenomenon of multiple unstable modes: if the unstable mode for k = 
has the squared frequency — \uj-\ 2 < then the k ^ modes will have an 
eigenmode with the eigenvalue k 2 /a 2 — |co>_| 2 which is positive for small a but 
will become negative whenever a > at = k/\u-\. The interpretation of these 
modes is unclear, there is no trace of them in the numerical simulations. 

With increasing a one expects the system to become unstable with respect 
to inhomogeneous, local bounces. This instability may be closely related to 
the sliding phenomenon for which a strong excitation of k ^ modes is 
crucial: the k ^ modes are of course spatially inhomogeneous; if they are 
excited substantially they may be interpreted as classical fluctuations. For 
k = n they have the general form of a chain of n bubbles, much in the way 
as they are depicted in Refs. [3J E] for the decay of the false vacuum at finite 
temperature. It could be interesting to investigate this relation between inho- 
mogeneous and homogeneous tunneling in a more quantitative way. Clearly, 
in view of this instability with respect to local bounces, our present analysis 
becomes unreliable at values a » 1. 
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The restriction to one compact space dimension was done in order to re- 
duce technical complications to a minimum. Similar phenomena are expected 
for compact three-dimensional spaces. The discrete spectrum of momenta 
then becomes a discrete spectrum of angular momenta, whose excitation 
will be small for small a and become more effective with increasing a. The 
potential seen by the homogeneous mode is again a double-well potential, 
so again resonances are expected to dominate the low a regime. If gravity 
is included, like for transitions in de Sitter space, renormalization becomes 
problematic, especially if quantum backreaction is included. Still it would 
be interesting to investigate along these lines the relation between homoge- 
neous transitions like the Hawking-Moss instanton and inhomogeneous ones 
like the Coleman-deLuccia bounce. 

Of course it would be even more interesting to follow, in real time, the 
local creation of real vacuum bubbles including the associated excitation of 
quantum modes with back-reaction. This seems to be out of scope with the 
presently available numerical methods and computer facilities. Our calcula- 
tions may elucidate in a qualitative way effects that are missed in the WKB 
approach to local transitions, even when quantum corrections are taken into 
account [111 ESI E3 • 
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